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Abstract: We provide a detailed analysis of the dynamics of moduli fields in the KKLT 
scenario coupled to a Polonyi field, which plays the role of a hidden matter sector field. It 
was previously shown that such matter fields can uplift AdS vacua to Minkowski or de Sitter 
vacua. Additionally, we take a background fluid into account (which can be either matter 
or radiation), which aids moduli stabilisation. Our analysis shows that the presence of the 
matter field further aids stabilisation, due to a new scaling regime. We study the system 
both analytically and numerically. 
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1. Introduction 

Theories beyond the standard model, such as string theory, predict the existence of massless 
(or nearly massless) scalar fields. These so-called moduli fields may contain, for example, 
the information about the dynamics of extra spatial dimensions, or, as it is the case of 
the string theory dilaton, they determine the coupling strength of gauge fields. In addi- 
tion, these fields usually interact, with gravitational strength, to other particles, thereby 
mediating a new ("fifth") force. Since so far we have neither observed any significant 
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time-dependence of gauge couplings nor detected new forces, the moduli fields have to be 
stabilised [1]. This should happen at some stage during the cosmological evolution, which 
is one of the problems addressed in string cosmology. It was emphasized in [2] that sta- 
bilising the string theory dilaton in a cosmological framework is a rather difficult, because 
of the the small barrier height separating a local Minkowski (or de Sitter) vacuum and 
the steepness of the potential. The situation can be improved if background fluids (either 
radiation or matter) or temperature dependent corrections are taken into account (for work 
in this direction see e.g. [3, 4, 5, 6]). 

In order to stabilise the moduli fields, a mechanism is required which gives them a 
mass. Additionally, the resulting vacuum energy should be very small and non-negative 
and the vacuum itself quasi-stable (i.e. with a life-time much longer than the age of the 
universe). Recently, compactification mechanisms with fluxes on internal manifolds have 
been extensively studied (see [7], a review can be found in [8]); a well-investigated scenario 
is the KKLT proposal [9]. Within this latter setup, moduli fields are stabilised and the 
resulting cosmological constant can be fine-tuned to be either zero or very small. To obtain 
a realistic scenario, a crucial ingredient in these mechanisms is a "lifting" procedure, which 
raised a supersymmetric AdS-vacuum to become a Minkowski or even a (quasi) de-Sitter 
vacuum. In the original KKLT proposal, an anti-D3 brane was added, explicitly breaking 
supersymmetry. In an alternative proposal, D-terms were considered to up-lift the AdS- 
vacuum, which requires the existence of charged matter fields [10, 11, 12]. 

Recently it was pointed out that interactions of the moduli fields with a hidden matter 
sector can also result in an effective uplifting to Minkowski or de Sitter vacua. The idea 
of this procedure is to take matter fields into account (such fields will be present in any 
realistic theory). The resulting interaction terms in the scalar potential lead to spontaneous 
supersymmetry breaking, with local minima which have a small and positive vacuum energy 
density This procedure is known as F-term uplifting (for recent work, see see e.g. [13, 14, 
15, 16, 17, 18, 19, 20] and references therein). An appealing feature is that this procedure 
can result in a small gravitino mass, while maintaining a high potential barrier, resulting 
in long-lived de-Sitter vacua [14, 17, 18, 21]. 

In this paper, we investigate the cosmological dynamics of the moduli fields and matter 
fields in the the context of F-term uplifing. As a toy model, we will couple the KKLT model 
with the Polonyi model [22], following [14, 20, 21, 23]. This model is simple enough to study 
the dynamical consequences of the interactions between the moduli fields and matter fields, 
but we believe that our results can be generalised to more complex models. In particular 
we are interested in whether the interaction will facilitate or impede the stabilisation of the 
moduli fields in the cosmological context. The paper is organized as follows: In Section 2 
we formalise the set-up and write down the essential equations. In Section 3, we discuss 
the properties of the potential. We review past results, in the context of our potential, 
in Section 4 and generalise to multi-field scenarios in Section 5 (two real fields). The full 
dynamics for two complex fields is presented in Section 6. We briefly consider the effects 
of a radiation background fluid in Section 7. We conclude in Section 8. We give further 
useful formulae in the Appendices and derive the new scaling regime. 
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2. Background Equations 



We begin by stating the four-dimensional N = 1 SUGRA action, which is of the form 

S = -JV^ (^-^R + Kfjd^^'d''^'^ + d^x, (2.1) 

where Kf^ — Q^i^^-j is the Kahler metric, are complex chiral superfields and V{<b) is 
the scalar potential, /^p is the 4-dimensional Newton constant, 

i4> = SttGn = 1. (2.2) 

The effective scalar potential is given, with given Kahler metric K and superpotential VF, 
as 

V = e^ (^K'~WiWD]W - 3WW^ , (2.3) 

where K'^^ is the inverse Kahler metric and DiW = diW + §§^W . 

The following equations of motion for the real and imaginary parts of superfields 
are [26] 

+ 3H^^ + ri.i^^^'k - ^j^'i) + ^K'^dj^V = 0, 

2 _ (2.4) 

^3} + 3H^} + r5,(v3i^| + ^'^^'j) + -K^'d,,V = 0, 

where (p^j^{(f\) refers to the real (imaginary) part of the scalar fields and dj^ (djj) are used 
to denote the derivative of the potential with respect to the real (imaginary) parts of the 
fields, respectively. The connections on the Kahler manifold are given by 

pn _ j^nl"^ (2.5) 

There are two different sectors in the model considered here; one involves a modulus 
field r, whereas the other sector contains a matter field C. Both fields will be taken as 
complex fields. The corresponding Kahler potential [21], which arises in type IIB and 
heterotic string theory, is 

= -31n(r + T) + |C|2, ^ ^ 

2.6) 

W =WiT)+W{C), 

Following [21, 23], for example, we consider the combination of the KKLT [9] and the 
Polonyi model, for which the superpotential is given by 



W{T) = Wo + Ae"^^ 
W{C) =c + fi'^C, 

so that the resulting scalar potential is given by 



(2-7) 
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V = 



3(r + Tf 



-Aae-^^iT + f ) - 3(Wo + ^e""^ + c + fi'^C)\^ 



+ 



H'C + {Wo + Ae-"^ + C + n'C) 



-3 



Wo + Ae-"^ + c + n'^C 



(2i 



The differential equations describing the dynamics of the field T are given by 

1.9.9 2T^ 
Tr + 3HTr - -{Tr - ) + -^Ot^V =0, 

2 . . 2T^ 
+ 3i7Ti - -TrTi + ^^T.F =0, 



(2.9) 



where we write T = Tr + iTi. Furthermore, we consider a background fluid with density 
and with equation of state 7 — 1 =^5/^5. Energy conservation dictates that 



+ 3H^pb = 0. 



(2.10) 



We consider C to be a complex field, C = Cr + iCi^ so that the equations of motion can 
be written as 



a + 3i/a + -5^1^=0, 

Q + 3TO + ^5c,y =0. 



The Friedmann equation reads 



(2.11) 



3^2 = ^(r2 + r, ) + {Cr + Ci) + F + 

Since is not a canonical field, it is convenient to define a new field 



(2.12) 



\nTr. 



(2.13) 



Additionally, the equations for the fields are easier to solve if we define a new set of variables 
as follows 



V 



X 



z = -e 



1 -.fl^Ti 



and 



p 



V3H' 

With these variables, the Friedmann equation becomes 

1 = + y2 + ^2 + p2 + + Qb, 



(2.14) 



(2.15) 



(2.16) 
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where fif, = pi,/^H^. The evolution of Hubble parameter H is 

H' = -^H [2x^ + + 2p^ + 2g2 + 7(1 - x"^ - y"^ - z"^ - - g^)] , (2.17) 
where the prime denotes differentiation with repeat to e-fold number N ^ defined by 

N = J Hdt. (2.18) 
The equations of motion for the fields are now given by 

x' = -3x + A^l^^ - 2^2 + [2x2 + 2^2 + 2p'^ + 2q^ + 7(1 - x'^ - - - - q^)] , 

+ 7(i-2;2-y'-^'-p'-9')] , (2.19) 

z' = -3z + ??y|y^ + 2xz + [2a;2 + 2z'^ + 2^^ + 2q'^ + 7(1 - x^ - - - - g^)] , 
and 



(2.20) 



p' = -3p + S^y'^ + ^p [2x'^ + 2^2 + + 2q'^ + 7(1 - x^ - _ ^2 _ ^2 _ ^2)j ^ 
= -3g + e^y^ + [2a;2 + 2z'^ + 2p'^ + 2q^ + 7(1 - x^ - y"^ - z^ - p^ - q^)] . 
In these equations, we have defined A and t] to be 

A = -l^ , = _,/?e^^l^ (2.21) 

For the remainder of this paper, we solve the above equations, in order to investigate the 
stabihsation of the moduh fields, and T^, in the presence of the matter fields, Cr and 
and the background fluid, p^. 

3. Properties of the Potential 

As a specific example in this paper, following [21], we take the potential parameters to be 

A=l, a =12, /i=10~^ c = /i^ (^2- a/s) . (3.1) 

Additionally, Wq is assumed to be real and is fine-tuned such that the potential is zero 
(or the cosmological constant A) at the local minimum. For the values above, the local 
minimum is found at the co-ordinates 

T^^^) - 3.3474 (0(1^) - 1.4797), T^^"^^ = 0, C^^^) - -0.7131, cf"^^ = 0. (3.2) 
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Figure 1: Scalar potential of the KKLT + Polonyi model in units of 10~^^tv~^. We plot the 
potential in the Tr-Cr plane (left hand side) and T^-T^ plane (right hand side), with the other fields 
set at their local minimum values shown in the text. 




Figure 2: Scalar potential of the KKLT + Polonyi model in units of 10~^^/^~^. We plot the 
potential in the (j)-Ti plane (left hand side) and (j)-Ci plane (right hand side), with the other fields 
set at their local minimum values shown in the text. 

In Figs. 1 and 2 we show the shape of this potential with the parameters given above. 
Since there are four independent fields, we fix two fields at their local minimum values and 
plot the remaining two fields. For convenience, we plot the canonical field, 0, instead of 
Tr- In Fig. 1 the local minimum of this model is seen, where (f) is stabilised at ^ 1.47. 
As ^ oc, the potential vanishes. Along the Cr direction, the potential increases as 
departs from the minimum. As can be seen, the minimum of Cr depends directly on the 
value of (/). Around the local minimum, Cr ^ —0.71. For (j) < —5, the potential can be 
approximated by 

V ^ ±±r^ (3.3) 



with Ti — Ci — where the minimum of is ^ 0. Note that the Cr direction is much 



-6- 





steeper than along the (j) direction. On the left hand side of Fig. 2, we show the sinusoidal 
shape of T^, with periodicity of 27r/a, with one minimum at zero. On the right hand side, 
the minimum of Ci is seen to exist at zero and shows symmetry about Ci ^ —Ci because 
we assume = in the figure. The values of A and 5 (defined in Eqns. (2.21-2.22)) over 
useful ranges are plotted in Fig. 3. 

It should be noted that, as the Ti and Ci evolve away from their minima, the shape 
of the potential changes from that shown in the figures. This will become important when 
we consider the evolution of the imaginary parts of the fields, C and T. 



4. Review of Single Field Dynamics 

We begin by reviewing known results [4, 25] and consider only real fields; we put the 
imaginary fields in the minimum point, = = 0, so that these fields do not evolve. 
This will enable isolation of the effects of each field (real and imaginary parts). Some of 
the results obtained will prove vital later. 

In this simplifed scenario, the potential becomes 



A^{3C^ + 4aTr{3 + aT^)} + Se^^'^^icCr + + ^^C^ + ujQ 

r ' 



24T3 

+ QAe"^^- {c{C'^ + 2aTr) + m^C^ + Cr{n^ + 2a/^^r^) + loC'^ + 2aioTr} 



(4.1) 



Further, to compare to single field dynamics, for negatively large 0, we set Cr in its local 
minimum {Cr — 0) initially, such that it does not evolve until the final stages of the 
evolution. Even with only (j) evolving, due to the background fluid and the potential, the 
dynamics are non trivial. 

In this section, we will closely follow [26] and will refer to an example trajectory 
shown in Fig. 4. We assume that (j) starts with zero kinetic energy from a point (j)f2 (this 
terminology will become useful later) and fi^ is close to unity initially (fi^o — 0.93). As the 
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5 10 15 20 25 30 35 40 5 10 15 20 25 30 35 

N 

(a) (b) 

Figure 4: Single Field Dynamics. Initial conditions are taken to be Q^o = 0.93, Cro = and 
(j)o = —15. The background fluid is matter, with 7 = 1. Note that Cr ip^) is zero until the last few 
e-folds. 



field starts to evolve, as long as A is small enough, the potential energy quickly dominates 
the dynamics, seen in Region 1 of Fig. 4(b). As the potential becomes more steep (A 
increases sharply around (j) ^ —5 for our potential), the evolution becomes kinetically 
dominated, called "kination" shown in Region 2. Once the background fluid dominates, 
the field is effectively frozen at a point (j)f^ (in Region 3), and the field only starts rolling 
again once this fluid has sufficiently diluted. 

As the field recommences to roll, one of two things can happen; either the scalar field 
dominates again (Region 1) or a scaling regime is found, in which the potential and kinetic 
energies track the background fluid (Region 4), which occurs for our trajectory. The choice 
of regime is specified by a scaling condition: 

37 < A^ (4.2) 

where we restate that 7 = 1 + ^5/^5. For our choice of potential, this condition is always 
satisfied when 0^—5 with a matter background fluid. The scaling regime ends when this 
condition is subsequently violated, which occurs if the potential flattens, for example near 
the minimum. If 0/3 < 0^^"^^ the end of scaling occurs near the minimum and the field 
stabilises. If 0/3 > (j)^^^^ but smaller than the local maximum value, the field will still 
stabilise. Conversely, if 0/3 is larger than the local maximum value, the field will roll to 
infinity. In our case, the fields are stabilised. 

In the following, we state the analytical solutions for the single field scenario. In 
particular, two initial branches can be distinguished; fi^o ^ 1 ^ind fi^o ^ 1 (the latter 
being described heuristically above). It is convenient to consider the new variables, x, y 
etc. (defined earlier) for analysis of the evolution. 
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4.1 ilbo < 1 

When flijQ <^ 1, following [24], the initial evolution is given by 



x' = xJ^y\ y'^-xJ^xy (4.3) 



where we assume A' ?^ 0, for which the solution is 



X 



= V(l-^^w)tanh (^A^livj , (4.4) 




y = V(l - Ubo) sech xJ -N . (4.5) 



This solution is valid until A' is no longer insignificant (0 ^ —5), at which point, the kinetic 
energy starts to dominate, 0.5 and kination starts. Due to the non-constancy of 

A, no approximation can be found for this period. However, during this stage, due to the 
insignificance of fi^, x^ increases to almost unity, as the kinetic energy fully dominates the 
evolution. Due to the expansion of the universe, the kinetic energy is diluted and slowly the 
background fluid becomes significant (il^ ^ 1, a::^ ^ 0), at which point the field becomes 
frozen at 0^3. 

When the potential is insignificant, the system during kination is greatly simplified 
and can be described by [4]: 

= -Sx + ^x [2x'^ +7(1- x^)] . (4.6) 
If X starts at a value xq, the solution is 



X, 







^ = ( 1 + L^e3(2-7)A^ ) (4.7) 



During this decay, the field moves a distance, A(/), given by 



^4> = ^n^M\^]- (4-8) 



3(2-7) 

after which the field freezes at <pf^. (Note that xq is never exactly unity since 7^ 0). 
4.2 fifto < 1 

When fi^o ^ 1? the the kinetic energy of never fully dominates, in that x'^ < 1^ before 
the background fluid takes over. As seen in Fig. 4(a), the background fluid is never totally 
negligible, which complicates the evolution. However, kination starts when A starts to 
increase { — 5, x^ ^ ^ 0.5). Once is negligible, moves a further distance given 
by Eqn (4.8) before freezing at . Although this regime cannot be solved analytically, 
one point can be made: in this case, the field freezes earlier (due to less kinetic energy 
overall, xq < 1), so that 0/3 is smaller than that when ^^^(O) <C 1. 

Numerically, we find that, for 7 = 1 and fi^o — 0.93 (an example that will prove useful 
later), stabilisation only occurs when (/)f2 ^ —17 and 0^3 < 1. We define these points as 

(f^min and (hm nx ' 
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5. Two Real Fields with a Matter Background Fluid 

This section extends the previous single field scenario, by considering the evolution of an 
additional scalar field. Specifically, we consider the evolution of (7^, in which direction the 
potential is much steeper than the (j)-fi.e\d previously considered. We will show that this 
leads to a new mechanism for stabilisation, due to the presence of a new scaling regime. 

We will explain the dynamics using several representative trajectories for which ^ 
shown in Fig. 5: 



I m = 


-23, 


Cr{0) 


= 10, 


Ob(0) 


= 0.01 


[solid line (blue)] 


II m = 


-15, 


CriO) 


= 3, 


Ofc(O) = 


= 0.01 


[dotted line (red)] 


III (p{0) = 


-23, 


CriO) 


= 3, 


nb{o) = 


= 0.01 


[dashed line (green)] 


IV (p{0) = 


-23, 


CriO) 


= 10, 




= 0.9 


[dotted line (magenta)] 



We take 7 = 1 (a background matter fluid) and consider two cases Jl^o = 0.01 and Jl^o = 0.9. 
The initial velocities of the fields are zero. 

Fig. 5 shows the evolution of the fields and Cr against e-folding number, A/", and the 
trajectories in the (/)-Cr plane. The evolution of each energy component for trajectories I 
and II are shown in Figs. 6 and 7. Trajectory III is identical to Trajectory II (although 
stabilisation does not occur) and is not shown. The evolution of the energy components 
for trajectory IV is shown in Figs. 8. Note that Fig. 6(b) strongly resembles Fig. 2 of [26]. 

As in the single field case, similar regimes can be identified: (1*) field domination, (2*) 
kination, (3*) freeze-out and (4*) scahng. The asterix (*) denotes two-field evolution, in 
order to separate it from the single field case. When Cr approaches zero, however, these 
regimes are followed by a period of oscillation around Cr = (Region 5*), further followed 
by the regimes (1-4) of single field evolution (as in the previous section). After this, the 
field is either stabihsed or rolls over. It should be noted that, due to the steepness in the 
Cr direction (5 > A), the Cr field arrives in its minimum quickly with the cj) field still 
running toward the minimum. Due to this, the second half of the dynamics follows single 
field evolution, as was discussed in the previous section. 

We will now derive some useful approximations for the first part of the evolution, when 
ICr-l 7^ 0. Depending on the initial value of Jl^, there are two branches of evolution, A and 
B. These are shown as a flowchart in Fig. 9, which also shows the path of each example 
trajectory. 

5.1 Branch A f^^o ^ 1 

Example trajectories for this branch are seen in Fig. 6 and 7, where il^o = 0.01. Initially, for 
a very brief time, the potential energy dominates but since the solution of field domination 
is not stable (see Appendix C) this energy is quickly transfered to the kinetic energy of 
both fields, (/) and Cr- This change of energy can be described by 




P 
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Figure 5: Various trajectories, each with differing initial conditions for (/), Cr and Q^o. Only real 
fields are considered (Tj = Q = 0). In (a), for Q^jq = 0.01, three trajectories are shown (I - solid 
blue, II - dotted red. III - dashed green). In (b), for l^^o = 0.9, only one trajectory is shown (IV 
- dotted magenta). The greyed region signifies initial conditions for trajectories which do not lead 
to stabilisation. The horizontal lines and regions i-iv are discussed in the text. In (c) and (d), the 
evolution of the fields are shown with respect to the e-folding number, for these four trajectories. 



for which the solutions are given by 



X = — ; tann 






p = — ; tann 

2/ = V(l - ^bo) sech 





(5.2) 



This is the extended solution of Eqns. (4.4) and (4.5) for two fields. The maxima of these 
solutions, which can be considered to define the start of kination, are approximately given 
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Figure 7: Breakdown of energy for Trajectory II with initial conditions 1^50 = 0.01, (j)o = —15, 
Cro = 3. We consider only the real parts of the fields (T^ = Ci = 0). 



by 



2 

^2 , 



A2 , 2 



2 ~0. 



y(2) 

Here we have neglected the evolution of fi^, but the error induced should be small. We shall 
use the notation that a subscripted bracket will denote the end of the regime of interest, 
here kination (Region 2*), but we omit the asterix for brevity. 
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(4*) 



;(5*); 



(1) ; (2) (3) (4) 

Traj IV 



10 20 30 40 50 60 70 80 



(b) 



Figure 8: Breakdown of energy for Trajectory IV with initial conditions 1^50 = O-Q, 
Cro = 10. We consider only the real parts of the fields (T^ = Ci = 0). 



= -23, 




Single Field 
Dynamics 



NO Stabilisation 



Stabilisation 



Figure 9: Flowchart of possibilities (Colors: Blue - Trajectory I, Red - Trajectory II, Green 
-Trajectory III, Magenta - Trajectory IV) 



During kination, we may assume 2/ — and the the solution is 

2 ^ ^(2) 

2^ 

(l-f^5(2))+f^5(2)e^(^--)(^-^(^))' 

' (l-f^5(2))+f^.(2)e^(^--)(^-^(^))' 
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When the potential energy neghgible compared to the kinetic energy, from the Klein- 
Gordon equation 

+ 3H(j) - 0, (5.5) 

for which the solution is oc a~^. Therefore the kinetic energy is proportional to a~^. This 
is the same for the Cr field. 

Alternatively, in this region, we can obtain 

X' X 



P_ _ P _ = constant, (5.6) 



where we have used the original definition of x and p for the third term. This gives a direct 
relation between x and p 

p = Cr - Cr^^^ = (0 - 0(2)), (5.7) 

where we have used the initial values of kination from Eqn. (5.3). The freeze-out point can 
be calculated by integrating these solutions 



= 



ACr = 



^1 - J^6(2) 3(2 - 7) 

P{2) 2V3 
^1 - 3(2 - 7) 



^6(2) / \ V ^6(2) 



(5.8) 



where A(f) = 0(3) — 0(2) and ACr = C'r(3) ~ ^^^(2) Note however, that when the initial 
background density is small enough, the transfer of potential energy to kinetic energy is 
so quick that we can consider starting at a point 00, C'rO with kinetic energies given by 

^(2),P(2). 

5.2 Branch B f^^o < 1 

An typical trajectory of this branch is shown in Fig. 8. If Jl^ is initially large, then the 
potential energy can never dominate and kination is not achieved. Instead, the equation 
of motion for the background fluid is given by 

n't, = [(7 - 2) (1 - Qb) + 2t/2] . (5.9) 

As potential energy is converted to kinetic energy, becomes negligible in a few e-folds 
and we quickly obtain the solution 

n,= (l + i^f^e^(7-2)ArA (5 10) 

Though the fields obtain some kinetic energy from the potential energy, it is suppressed 
compared to the background density. Quickly, fi^ ^ 1 and the fields become frozen. This 
process is very fast and we can make the approximation 0(3) ^ 0o and C^(3) ^ Cro 
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5.3 Scaling 

We calculate the full scaling solution and condition in Appendix B. We will use the main 
results in this section. The two-field scaling condition is given by 

37< A^ + 5V2. (5.11) 

At the freeze-out point, if the condition for scaling is satisfied, then the fields track the 
background fluid. We define 

r«^^. (5.12) 

Making the assumption that 

- Tc.a - r^cr - 1, (5.13) 
which is reasonable in the ranges 4> < —8 and IC^-I <; 1 for our potential, we obtain 



Xn — 



where 



7=7(l + 5^j . (5.15) 

This is a valid solution only when Tij ^ 1. When (7^ = 1, however, 6 changes form and 
TcrCr begins to decrease (and later increases again). This deviation from F ^ 1 leads to 
an increase in fi^ (as seen around ^ 35 in Fig. 6). Scaling proceeds until Cr goes to 
around 0, where the scaling condition is violated. 

During scaling, the background fluid density takes the form 

n, = l-xl-pl-y'', = 1-^ = 1- j^^^ (5.16) 

and the trajectory is specified by 

- = ^ = (5-17) 
Xc V2X d(f) 

Prom Fig. 3, we may take A approximately constant and from Eqn. (4.1) 

5 « -2{Cr + ^) (5.18) 



SO that 



^ ^ + W(3) 
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At the end of the scahng regime, as found numericaUy, the fractional energy of a matter 
background fluid is 

- 0.93, (5.20) 

and the kinetic energy of (f) is small. We find that this fractional background density is 
numerically independent of the initial values of the dynamics and is determined only by 
the existence of scaling solution (i.e. potential dependent). Additionally, from Eqn. (5.19), 
when the trajectory reaches = 0, we find 

0(5) ^ 0(3) - ^ log . L . (5.21) 

where either [0(3),C^(3)] are calculated from Eqn. (5.8) (Branch A) or [(t)o,Cro] (Branch B). 
For a matter background and fi^o ^ 1? trajectories always scale. When fi^o ^ 1? oiily 
trajectories with IC^-ol ^ 4 have time to enter a scaling regime. This region is depicted by 
regions (i) in Fig. 5(a) and separated with a dotted line. 

5.4 Oscillation 

If the trajectory undergoes scaling, then it is clear from Eqn. (5.14) that, for our potential 
at least, Pc > Xc and there is more kinetic energy in Cr than in (j). When the trajectory 
reaches the Cr minimum {Cr = 0), the trajectory oscillates but is almost constant; it 
acts as if it is almost frozen at a point (f)f2 with ^ 0.93. 

If the trajectory does not scale (like Trajectories II and III), the trajectory gently 
oscillates around the minimum, with not fixed. Due to the relative gradients, A and 5, 
the kinetic energy is shared between the fields, but the expansion of the universe acts to 
damp the oscillations. Finally, the background fluid dominates and the fields stop rolling 
at a point 0^2- At this point, we find numerically that fi^ ^ 0.9 once again. 

After freeze-out, Cr ^ and the trajectory enters single field dynamics as discussed in 
Section 4. 

5.5 Stabilisation Regions 

We shall characterise a stabilising trajectory by the initial conditions for and Cr which 
lead to final stabilisation. Samples of the stabilisation regions for two real fields can be seen 
in Fig. 5. Greyed out regions denote initial field values for trajectories which do not result 
in stabilisation. We will see that the important factor in stabilisation is the point 0^2, the 
value at which the 0- field freezes on the Cr = line, after either scaling or oscillation. 

5.5.1 n^o < 1 

We will first explain the stabilisation region for a large background fluid density, such as 
when fi^o = 0.9 as shown in Fig. 5(b). A trajectory starting at large \Cr\ will follow similar 
dynamics as Trajectory IV. After freezing at a point [0/i,C'^(/i)], the fields will quickly 
reach a scaling regime. When Cr ^ 0, scaling ends and the Cr^-field oscillates around the 
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"valley" and (j) is effectively fixed at 0j2- Henceforth, the trajectory follows single field 
dynamics, with initial conditions (/)f2 and 0.93. 

As discussed in Section 4, stabilisation depends merely on these initial conditions. 
Stabilisation may occur only when —17 < 0/2 < 1 for il^ = 0.93, as is seen in Fig. 5(b). 
The slight asymmetry of the stabilisation region (seen for large 0) is due to the asymmetry 
of the potential; the local minimum does not occur at Cr = 0. Trajectories which start 
close to the minimum are easier to stabilise. 

We have numerically simulated trajectories only within \Cr\ < 10. Note however, that 
all trajectories with within the stabilisation bounds at C^- = will lead to stabilisation. 
Therefore the allowed Cr region extends to \Cr \ oc, due to the scaling solution. 

5.5.2 Qto < 1 

The stabilisation region for fi^g <^ 1 looks vastly different from the previous example, as 
can be seen in Fig. 5(b). However, the underlying structure is the same. In regions (i), 
similar constraints are seen. These are due to scaling solutions, which result in the (/)-field 
freezing on the Cr axis between — 17 < 0/2 < 1 as before. 

In region (ii), however, though the scaling condition is satisfied, the field does not have 
enough time to reach this scaling regime before reaching Cr = 0. Thus the field has enough 
kinetic energy at the crossing point to oscillate around = (seen in Fig. 5(b)). As 
described earlier, the trajectory settles into the minimum at = at a point 0/2- As 
long as (f)f2 is within the range — 17 < 0/2 < 1, the trajectory will lead to stabilisation. 
Therefore Trajectory II stabilises (0/2 ^ —12) whilst Trajectory III does not (0/2 !^ —20). 

The conical boundary of region (ii) depends on the energy density of the background 
fluid. For smaller fi^o? the kinetic energies are less damped and oscillate more widely. Since 
the field gains more kinetic energy initially (from potential energy) and is less damped, the 
fields travel further before settling at a point 0/2. Therefore as Q50 decreases, the conical 
bound moves left and the stabilisation area increases. 

The asymmetry seen between regions (iii) and (iv) is once more again to the asymmetry 
of the potential; trajectories starting in region (iii) are pushed towards the minimum (due 
to the gradient) and stabilise, whereas trajectories from region (iv) are pushed away from 
the minimum. 

6. Dynamics of Complex Fields 

In this section, we include the dynamics of the imaginary parts of the fields, Ci and T^, in 
order to fully generalise our model. The existence of non zero imaginary parts changes the 
shape of the potential in the (f)-Cr plane during evolution. As an example we consider a 
trajectory with initial values: 

= -2, a = 0.25, Ti = 0, a = l, f^60 = 0.9. (6.1) 

The evolution of the fields for this trajectory are shown in Fig. 10. It can be seen that, 
even though T^o = 0, this field is given a "kick" since Ci ^ We show the evolution of each 
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Figure 11: The evolution of each energy component for the trajectory given by the initial conditions 
in Eqn. (6.1). (a) The breakdown of kinetic energies, (b) The total breakdown. Note that the 
kinetic energy of Ti is negligible since T^o = 0; it is not a generic result. 

energy component in Fig. 11, where the kinetic energy of Ti is seen to be negligible. This 
is due to the fact that T^o — and is not generic. 

In Fig. 12 we show the contour plots of the potential in the (j) — Cr plane for the above 
trajectory at e-folding numbers N — 26, 27, 29, 30. In the top lefthand plot, where Ti :^ 0.1 
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N=26 Tj= 0.131 Cj= 0.703 N= 27 Tj= 0.110 Cj= 0.665 




Figure 12: The potential at different efolding values {N = 26,27,29,30), for two complex fields 
{Ti and Ci are non-zero), showing how the existence of the minimum (in the (j)-Cr plane) depends 
on the imaginary field values. The initial conditions are given in Eqn. (6.1). At N ^ 26, the local 
minimum has almost disappeared since Ci is large, but reappears at later times when Q ~ 0. The 
evolution of the potential works to nudge the fields into the local minimum. 

and Ci ^ 0.7, the local minimum almost disappears. However, at larger N and smaller C^, 
the potential changes shape and finally around ^ the local minimum re-appears and 
the fields find this minimum. Thus we see that evolution of the imaginary fields can change 
the possibility of stabilisation. In fact, the evolution of the potential works to push the 
fields into the local minimum; as the fields get close to the position of the local minimum, 
the potential barrier effectively increases (the plateau is raised) and the fields are pushed 
back towards the re-emerging minimum. 
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Table 1: Table showing the stabihsation measure, Aah^ (percentage of the total area which leads 
to stabilisation) as shown graphically in Figure 13. These numbers are indicative of trends only, as 
explained in the text. 

In order to understand the full stabilisation region, a coarse 4-dimesional grid of tra- 
jectories was analysed numerically, in which each point represents the initial conditions for 
[(j)^Ti^Cr^Ci]. Three such grids were run, for = 10-1°, 10-2, 0.9. The results are shown 
in Fig. 13 for the 3-D volume (j)-Cr-Ci. We do not show the direction, since all values 
lead to identical surface plots (this is due to the symmetry seen in Fig. 2). 

It is convenient to define a measure of goodness of stability. Firstly, to determine 
whether the background fluid adds stabilisation, we define an overall normalised, 4-D vol- 
ume, V, which is given by the volume of stabilised trajectories as a percentage of the total 
volume. In addition, we may calculate the ratio of stabilised trajectories to the total tra- 
jectories on a given 2-D plane (again, as a percentage), Aab^ where a and h denote the fixed 
values of the remaining fields. (For example Aci=i{),Ti={) represents the stabilisation area 
of the (t)-Cr plane when = 10 and T^o = 0). 

The numerical results for the measure of goodness are given in Tables 1 and 2. It 
should be well noted that the numbers quoted are intended to show a trend and are only 
indicative. Numerically, we rely on a finite volume, with finite (coarse) spacing. We chose 
a grid size that was numerically feasible within timescales. The may preclude us from 
resolving all the fine structure. We do note, however, that the fine structure arises from 
the asymmetry of the potential and the fine-tuning of the trajectories. This will be more 
relevant when we consider a radiation background fluid. 

Our results indicate that a background matter fluid aids stabilisation, consistent with 
previous work [4, 25, 26]. In addition, the presence of matter fields further aids stabilisation. 
Importantly, while a non-zero value of Cr assists the stabilisation of 0, the imaginary part 
Ci further increases the stabilsation regions. 

7. Radiation (7 = |) 

For completeness, we briefly consider the effects of a radiation background fluid, in place 
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Figure 13: Matter (left column) and Radiation (right column), from top to bottom: Q^q = 
10~^^, 10~^, 0.9. A measure of stabilisation is the percentage of the total area which leads to 
stabilisation. Values are given in Table 1. 

of a matter fluid. All the equations in the preceding sections are valid, where 7 = 4/3. 
The results are shown in the righthand column of Fig. 13 and Tables 1 and 2. The main 
conclusions are: 

1. For all scenarios, the stabilisation region is much reduced, compared to when 7=1. 
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Table 2: Table showing the stabilisation measure, V, (percentage of the total 4-D volume which 
leads to stabilisation). 2-D cross-sections are shown graphically in Figure 13. These numbers are 
indicative of trends only, as explained in the text. 

2. Increasing fi^o again aids stabilisation. 

3. Due to the reduced cosmological friction acting on the fields, more fine-tuning is 
necessary for trajectories to find the local minimum. 

4. We observe the same trend of stabilisation when Cr and Ci are non-zero as for the 
matter background fluid, although the effects are less pronounced. 

8. Conclusion 

In this paper, we have considered moduli stabilisation in the KKLT scenario coupled to a 
Polonyi matter field, a model which has been previously studied in the context of particle 
phenomenology. In this scenario, the AdS vacuum is uplifted by F-terms in supergravity. 
We have examined the evolution of two complex fields (one modulus and one matter field) , 
taking also a background fluid into account, either matter or radiation. We find that the 
presence of both the background fluid and the matter field enlarge the region leading to 
stabilisation of the moduli fields, due to a new scaling regime. Although a more detailed 
treatment is necessary, we believe that our conclusions will hold for similar scenarios, 
where the new direction is steep. A matter background fluid aids stabilisation more than 
radiation, consistent with previous studies. 
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Appendix 

A. Useful Formulae 

From the definition of A, 5, 77 and 0, one obtains 

J = VgxX [1 - r^^] + Vgztj [1 - p^tJ + V3p6 [1 - r^cr] + [i - r^cj , 

J = VGxX [1 - Cr] + V6zr/ [1 - Fc.tJ + VSpS [l - Tc^Cr] + VSqO [l - Fc.cJ , 
^ = 2:r + V6xX[l - F^tJ + Vezr][l - Ft.t,] + V3p6[l - Ft,c.] + V^qO[l - Fr.cJ, 
J = VGxX [1 - F^cJ + Vezrj [1 - Fc.tJ + V3p6 [l - Fc.cJ + VSqO [l - Fc.cJ , 
where we define 

We take a change of variables 

_ 1 _ 1 _ 1 _ 1 

and we define X, Z, P and Q to be 

x = e^X, y = e^Y, z = ct^Z, p = eCrP, Q^^dQ- 
With these new variables we obtain: 

+ 7(1- - _ ^2^^2 _ ^2^p2 _ ^2^^q2) ] 
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and 



X' = ^X-?,X+ \IIy^ - 2^Z2 - 

A V 2 e^fy H 

z' = iz-zz^f-[ltJy^^2.,xz-jLz 



Finally, from above, we find 

^ = - {^/6X [1 - r^^] +V6Z[1- T^T,] +V3P[1- F^Cr] + VSQ[l - T^cJ } (A.l) 
-[V6X[1- r^cr] +V6Z[1- Tc^T,] +VlP[l- Tc^cr] + V3Q[l - Tc^c, 



(A.2) 

^ = - |2e^X + V6X [1 - T^T,] +V6Z[1- Tr^rJ +V3P[l- Tr^Cr] +V3Q[1- Tt.cJ | 

(A.3) 

e'c, = -ec, { ^/6X [1 - r^cj +V6Z[1- T^tJ + a/SP [l - Tc.cJ + x/3Q [l - Tc.cJ } 

(A.4) 

B. Tracking Solution 

In the following, we generalise the scaling regime found in [24] for two real fields. We 
therefore consider = = in the following. The simplified equations of motion 
considered (with Q — Z — {)) are 

x'= -y6(r^^-i)x2 _^(r^c^ - i)xp - 3x + Y^r^ ^^ .^^ 

+lx [2(4X^ + elP^) + 7(1 - _ ^2^2 _ ^2^^p2)] 

Y' = -VQ{T^^-1)XY -V^{T^Cr-^)yP-\l\xY -^PY (B.2) 

+^Y [2(e|x2 + elP'') + 7(1 - e^X^ - e^yS _ ^2^^p2^] ^ 



.2 



p' = -V6(r^a - 1)^^ -V3(rc.c. - i)p' - 3P + ^^i^' (b.3) 



2 



+|P [2(4X2 + P2) + 7(1 - 4x2 - 4y2 _ ,2^^p2)] 
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Small and or Fij « 1 imply that is almost constant, as seen in Equations (A.l) and 
(A. 2). Therefore the "instant critical point" is obtained from solving X' = Y' = P' = 0. 
Prom Eqns. (B.1)-(B.3), we obtain 



(B.4) 



^'(^ - 754;^) = 2(r00 - ^4>Cr)X^P + ^(^<PCr - rCrCr)XP\ (B.5) 

It is reasonable to assume (in the regions we consider in the paper) 

^ ^CrCr — r^Cr — ^5 

from which we find 



-X. 



Plugging this solution into Eqn. (B.4), both P = P{X) and Y = Y{X). Solving Eqn. (B.l) 
(using X' = 0) leads to 



2A' 



Xr, = M--, Pc = ^^,, yc = 14 ( 2-7- 1 , 



2A2 



where 



7 = 71 + 



2A2 

C. Stability of the critical points 

We expand about the critical points 

X = Xc + u, Y = Yc + v, 
which yield, to first order, the equations of motion 



(B.6) 



(B.7) 



P = Pc + W, 









v' 


= M 


V 


\w') 




[wj 



Assuming F's are 1 (this assumption is quite good for the region of approximately < —8 
and \Cr\ > 1), we find the eigenvalues, m^, of M. We also assume < 7 < 2 for the 
baryotropic fluid. 
Fluid- dominated solution 

For this solution, Xc = Pc = Vc = 0- We find a saddle point for < 7 < 2: 

3 3 



mi = -7, 



m2 = rris 



:(2-7)- 
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Kinetic- dominated solutions 



Here, xl+ = 1 and yc = 0. We find an unstable node for Xxc + {S/y/2)pc < a/6 and a 
saddle point for Xxc + {6/^/2)pc > \/6: 

mi = 0, m2 = 3(2-7), ^ \/f ~ '^^^ ~ ^) * 

Scalar field dominated solution 

In this case, Xc = = Vc — ^ ~ T2 (^-^^ +5^)- There is a stable node for 

+ 52/2 < 37 and a saddle point for 87 < + < 6: 

mi,2 = -3 + ^{S^ + 2A'), ms = -87 + \{S^ + 2A'). 

Scaling solution 

For the scaling solution, Xc, and are given in Eqns. (B.6) and (B.7). We find a stable 
node for 87 < A^ + 5^/2 < 247^/(97 - 2) and a stable spiral for A^ > 247^/(97 - 2): 



mi 



^2,3 



-2(2-7), 



l± Wl- 



87(A2 + Sy2 - 37)) 



(A2 + 5V2)(2-7) J 
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